Electron screening and excitonic condensation in double-layer graphene systems 



Maxim Yu. Kharitonov 1 and Konstantin B. Efetov 1 ' 2 
Theoretische Physik III, Ruhr-Universitat Bochum, Germany 
2 L.D. Landau Institute for Theoretical Physics, Moscow, Russia 
(Dated: August 15, 2008) 

We theoretically investigate the possibility of excitonic condensation in a system of two graphene 
monolayers separated by an insulator, in which electrons and holes in the layers are induced by 
external gates. In contrast to the recent studies of this system, we take into account the screening 
of the interlayer Coulomb interaction by the carriers in the layers, and this drastically changes the 
result. Due to a large number of electron species in the system (two projections of spin, two valleys, 
and two layers) and to the suppression of backscattering in graphene, the maximum possible strength 
of the screened Coulomb interaction appears to be quite small making the weak-coupling treatment 
applicable. We calculate the mean-field transition temperature for a clean system and demonstrate 
that its highest possible value T"™ ax ^ 10~ 7 ef < 1 mK is extremely small (ep is the Fermi energy). In 
addition, any sufficiently short-range disorder with the scattering time r < fi/T™ ax would suppress 
the condensate completely. Our findings renders experimental observation of excitonic condensation 
in the above setup improbable even at very low temperatures. 

PACS numbers: 73.63.-b, 72.15.Rn, 81.05.Uw 



INTRODUCTION AND MAIN RESULT 

The possibility of excitonic condensation (EC) in 
metallic systems was originally proposed by Keldysh and 
Kopaev [lj for semimetals with overlapping conduction 
and valence bands. They have shown that the attractive 
Coulomb interaction between electrons and holes leads to 
an instability towards formation of bound electron-hole 
pairs analogous to the Cooper instability in supercon- 
ductors. Somewhat later, it was suggested 0] that EC 
could be realized in a double-layer system of spatially 
separated electrons and holes. Experimental efforts to- 
wards the observation of EC were mainly concentrated on 
semiconductor double quantum well systems 0, 0, [!, 0] 
and experimental data speak in favor of the existence of 
EC in electron-hole bilayers 0,0, 0] and electron-electron 
bilayers in the quantum Hall regime 0] . 

Since the carrier density in graphene, including its po- 
larity, can effectively be controlled by various means, 
graphene-based systems may also seem attractive for the 
realization of EC. Indeed, several ideas on how one could 
obtain EC in graphene have been suggested recently. One 
possible way to create interacting electrons and holes is 
to apply a strong in-plane magnetic field to a single layer 
of graphene [7J ■ Such a magnetic field acts on the spins of 
the carriers only, and the Zeeman splitting creates elec- 
trons with one spin polarization and holes with the op- 
posite polarization in an initially neutral sample. A de- 
tailed theory of EC in such a setup has been developed 
in Ref. 0]. 

A double-layer graphene system (Fig. [T]) as a candi- 
date for the observation of EC was proposed recently in 
Refs. 0, @, [nj. If two graphene layers are separated by 
an insulator, electrons in one layer and holes in the other 
can be obtained by applying external gate voltage. Rel- 
atively high values of the Fermi energy cf ~ 0.3 eV that 
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FIG. 1: Excitonic condensate in a system of two spatially 
separated graphene layers. Electrons and holes in the layers 
are induced by applying the external gate voltage. 
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can be achieved in graphene by using gates 
obvious advantage, since ep serves as a high-energy scale 
of the effect in such a setup. Solving the gap equations 
numerically, the authors of Refs. 0,11(1 provided an es- 
timate Tbkt ~ O.Icf for the critical temperature of the 
Berezinski-Kosterlitz-Thouless (BKT) transition and ar- 
gued that Tbkt could thus reach room temperatures. 

However, in the analysis of Refs. 0, [nj, the screen- 
ing of the Coulomb interactions by the carriers in the 
graphene layers was not taken into account. Clearly, the 
two-dimensional screening cloud formed around a probe 
charge in a graphene sheet screens the field of the charge 
both in the off-plane and in- plane directions, although 
not identically. Therefore, screening affects not only the 
intralayer but also the interlayer Coulomb interaction in 
the double-layer setup. 

In this paper, we demonstrate that taking screening 
into account is essential as it drastically reduces the tran- 
sition temperature compared to the estimate obtained 
in Refs 0, [1(3] neglecting screening. In fact, screening 
sets the upper bound for the interaction strength, yield- 
ing for the maximum possible value of the dimensionless 
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Here N = N s N v Ni = 2 3 = 8 is the total number of 
electron species in the system originating from two pro- 
jections of spin (N s = 2), two valleys (N v = 2), and two 
layers (iV; = 2). Moreover, the chiral nature of quasipar- 
ticles in graphene leads to the suppression of backscat- 
tering. Consequently, the maximum interaction strength 
A™ that determines the transition temperature appears 
to be actually two times smaller than A max , 
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As follows from Eqs. ([T]) and ([2|), the large num- 
ber of electron species and suppression of backscattering 
make the maximum possible value A™ ax of the interac- 
tion strength numerically quite small. This justifies the 
applicability of the weak-coupling BCS approach to the 
problem, since 1/16 can safely be considered as a small 
parameter. As a result, for the highest possible value of 
the mean-field transition temperature we obtain 
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exp(-l/A™ ax ) e F = exp(-16) e F 
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e F . 



(3) 

This is the highest possible value of the critical tempera- 
ture of the excitonic condensation that could be achieved 
in a perfectly clean double- layer graphene system. Enter- 
ing Eq. ([3]) as the exponent, the large value of the inverse 
interaction strength 1/A™ ax = 16 results in a drastic re- 
duction of the transition temperature. In order to achieve 
the maximum value ((3]), the interlayer distance d must 
be much smaller than the Debye screening length x^ 1 , 

2nd < 1. (4) 

The most optimistic estimate would thus be T™ ax ~ 
1 mK for ep ~ 0.3 eV and would require d < 0.2 nm. 



CALCULATIONS 

We now present the details of derivation of Eq. 
The bare strength of Coulomb interactions in graphene 
is not that small. For Si02 used as an insulator embed- 
ding graphene sheets, typical values of the dimensionless 
coupling constant are r s — e 2 / (ev) ~ 1 (e is the dielectric 
constant of the insulator and v is the velocity of the Dirac 
spectrum, we set h = 1 throughout this section). This 
questions the applicability of the weak-coupling approach 
to the problem of EC, suggesting, at the same time, that 
the transition temperature could be quite high 

It is known, however, that in a fermionic system with a 
large number N 3> 1 of independent fermionic species the 
interactions are significantly weakened. Physically, large 



N makes screening very effective, since all N species par- 
ticipate in the screening of interactions between fermions 
of each particular species. Screening reduces the coupling 
constant from r s to the value X/N < 1: r s — > 1/N. Ef- 
fectively, the system becomes weakly interacting, despite 
of the fact that the bare Coulomb interactions may be 
not weak (r s > 1). 

The large-iV approximation was already used for a 
single-layer graphene @, E2, E3] before. In a single layer, 
the number of species is equal to Ni = N S N V = 4 due 
to two projections of spin (N s = 2) and two valleys 
(N v = 2). This value is not exceptionally large, but does 
give hope that the large- N approach adequately describes 
graphene physics. In a double-layer system the situation 
is better: since each electron can belong to either one 
of the layers, one has an additional "which-layer" degree 
of freedom (Ni = 2) making the total number of species 
N = N s N v Ni = 8. It would already be quite reasonable 
to treat N — 8 as a large parameter. Therefore, large- 
TV approximation seems to be particularly suitable for a 
double-layer graphene system and is expected to provide 
good quantitative predictions. 

Below we employ the large- N approach to the double- 
layer graphene system (Fig. [1} treating TV = 8 as a large 
parameter and calculate the mean-field critical tempera- 
ture Tr of EC. The calculations follow closely those of 
Ref. [7|. Of course, the mean-field treatment of two- 
dimensional systems is not necessarily a good one due 
to strong thermal fluctuations. It is, however, sufficient 
for our purposes, as our main goal is to demonstrate that 
already the mean-field critical temperature is extremely 
low. The temperature Tbkt of the actual BKT transi- 
tion can only be lower than the mean-field T c we calculate 
here. 

Within the large- TV approximation, the diagrammatic 
series for the effective interaction between the electrons 
is identical to that of the random-phase approximation 
(RPA), which describes linear screening. For the prob- 
lem at hand, the relevant transfer momenta q are in the 
range < q < 2p F , with p F — e F /v the Fermi mo- 
mentum, and one has to use the exact expression for the 
polarization operator (and not its limit form for q <^ p F ). 
However, the static polarization operator in graphene 14 1 
does not depend on momentum at all in this range and 
equals H(ui = 0, q) = N s N v i/, where v = e F /(2irv 2 ) is 
the density of states per one valley and one spin. As a 
result, for the screened interlayer Coulomb interaction in 
the momentum space we obtain 



V(q) 



27re 2 exp(— qd) 



q + 2x + x 2 [l — exp(—qd)]/q 



q<2p F . (5) 



In Eq. ||SJ), q is the absolute value of the in-plane two- 
dimensional wave vector, d is the distance between the 
layers, e* is the effective electron charge screened by the 
insulator embedding graphene sheets, e\ = e 2 /e, and 
x = 2-kN s N v e\v is the inverse Debye screening length 
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FIG. 2: The screened V(q) [Eq. ([5}, solid line] and unscreened 
V UB {q) [Eq. (|8]), dashed line] interlayer Coulomb interaction, 
responsible for the excitonic instability. The values d = and 
r s = e 2 /(ev) = 1 were used, ^ is the density of states. At 
relevant momenta q ~ pf, the unscreened interaction poten- 
tial overestimates the actual screened one by about 10 times. 
The screened Coulomb interaction reaches its maximum at 
q — 0, the universal value ^ of which is achieved for 2xd 1 
[Eq. Q]. 



in each layer. We assume the same Fermi momenta pf 
of electrons and holes (this can be achieved by tuning the 
gate voltage), since any difference between them would 
be suppressing the condensate in a way Zeeman splitting 
suppresses s-wave superconductivity. 

The screened Coulomb interaction © is a decreasing 
function of q (Fig. [5]) and reaches its maximum at q = 0, 



V(q = 0) = 



2nei 



2x + x 2 d' 



(6) 



The maximum of Eq. ([5]) is achieved, if the distance 
d between the layers is smaller than the Debye ra- 
dius [Eq. Q], and equals 



V max (q = 0) 



1 



(7) 



The factor 2 that enters the denominator in Eq. (0 is 
due to "which layer" degree of freedom (JVj = 2), since 
each carrier can belong to either one of the layers. Equa- 
tion (J7|) leads to Eq. |T]) for the maximum value of the 
dimensionless coupling constant A max = vV m&x {q = 0). 

In contrast to the above calculation, in the analysis of 
Refs. 0, H3] the unscreened form 



V us (q) = 2irelexp(-qd)/q 



(8) 



of the Coulomb interaction V(q) [Eq. (J2J)] was used, see 
Eq. (4) in Ref. [j| and inline formulas before Eq. (2) in 
Ref. As seen from Eqs. © and © and Fig. [1 the 
unscreened form V us (q) is valid for q 2k = Nr s pp, 
but significantly overestimates the actual screened inter- 
action V(q) for relevant momenta q < 2pp. For the value 
r s = 1 used in Ref. [§| and typical for Si02 as an insu- 
lator, one obtains V us (Pf) ~ 9V(pf)- Using the un- 
screened form of the Coulomb interaction in Refs. 0, [l(| 



resulted in the estimate Tbkt ~ 0.1 ep and, as it ap- 
pears, led an overestimation oITbkt by a factor 10 6 , see 
Eq. ©. 

In order to obtain the mean-field transition temper- 
ature T c , we derive the linearized gap equation for the 
order parameter 



A(r-r') = V(r-r')(0 e (r)^(r')). 



(9) 



Here, V(r — r') is the interaction (|5|) in the coordinate 
space and e ,/i are the Dirac spinor fields of electrons and 
holes in the graphene sheets. The matrix structure of the 
order parameter in the sublattice space is predetermined 
by chirality, but can be arbitrary in the valley and spin 
spaces. Using the standard BCS approach, we arrive at 
the linearized gap equation 



A(n) = v\m 



££ 
T 



dn' 
~2T 



V(pF\n-n'\)V(n')A(n')V(-n'), 

(10) 

where the two-dimensional unit vectors n and n' repre- 
sent the direction of the electron momentum and 'P(n) = 
(1 + rn)/2, with t x and r y the Pauli matrices in the sub- 
lattice space. The value of temperature T, at which a 
nonzero solution A(n) to Eq. (fl"0|) appears, determines 
T c . Solving Eq. (10]), we obtain 



T c sa exp(-l/A c ) e F , 



where 
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2pp sin - 
2 
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(12) 



The exact numerical value ~ 1 of the prefactor in Eq. (fTT|) 
cannot be obtained within the logarithmic accuracy of 
the mean-field approach. The form of Eqs. (fTTjl and 
(fT2f and the matrix structure of the solution A(n) in 
the sublattice space are identical to those in Ref. Q [see 
Eqs. (5.23)-(5.31) therein]. At the same time, the form 
([5]) of the interaction V (q) is different here. 

The maximum possible value A™ ax [Eq. ©] of the in- 
teraction constant A c [Eq. (fT2|) ] and, thus, the highest 
possible transition temperature T™ ax , [Eq. ©J, are ob- 
tained by inserting Eq. ((7|) into Eq. (fTSj) . This corre- 
sponds to the limit p F <C 2k <C 1/d [Eq. @], where 
the condition pp <C 2k is automatically satisfied, since 
r s - 1 and 2k/ p F = r s N > 1. The factor (1 + cos6>)/2 
entering Eq. (|12[) is a consequence of chirality. It sup- 
presses backscattering and reduces A™ ax by a factor 2 
compared to A max [Eq. ©], see Eq. ©. 

The obtained small value of AJ. nax = 1/16 justifies 
the very applicability of the weak-coupling BCS ap- 
proach to determining T c , within which the logarithm 
ln(ep-/T) « 16 has to be large. Therefore, the criti- 
cal temperature T c does exponentially depend on the in- 
verse coupling constant 1 / A c , which leads to its extremely 
small value [Eq. ©]. 
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FIG. 3: Sensitivity of the excitonic condensate to the impu- 
rity scattering. Impurities with the size of potential smaller 
than the interlayer distance scatter electrons and holes not 
identically, thereby breaking electron-hole pairs and suppress- 
ing the condensate. 



DISCUSSION AND CONCLUSION 

Let us now discuss the obtained results. Remarkably 
enough, as Eqs. H])-© demonstrate, the specifics of the 
graphene spectrum (chirality and valley degrees of free- 
dom) appears to be very unfavorable for the realization 
of EC in graphene-based devices. At the same time, 
this is not so for double-layer systems based on materi- 
als with "conventional" metallic spectrum, such as, e.g., 
GaAs/Alj;Gai_a;As heterostructures used so far experi- 
mentally d, i, Hi]. 

Indeed, in such systems the max- 
imum interaction strength A' c nax = l/(N s Ni) — 1/4 is 
four times larger than in graphene [Eq. ([2])] due to the 
absence of the valley space and chirality. This value is not 
that small and the system could be on the verge of the 
weak-coupling limit. Therefore, one does not get such a 
small value of the transition temperature as we obtained 
for graphene. 

It is also instructive to mention for comparison that 
in a single-layer graphene subject to the in-plane mag- 
netic field 0] one obtains A™ ax = 1/(2 N S N V ) = 1/8 for 
the interaction constant, since the system consists of only 
one layer (JVj = 1) [see Eq. (5.29b) in Ref. Q], and the 
exponential factor exp(— A™ ax ) = exp(— 8) « 3 • 1CP 4 in 
Eq. ([3]) is not as small. However, the Zeeman splitting 
energy enters Eq. ([3]) instead of £f , which cannot be ex- 
tremely high even for experimentally very high magnetic 
fields B. For B « 40 T one can estimate T c max - 20mK. 

There is another factor that is unfavorable for the re- 
alization of EC in double-layer systems. Namely, the 
excitonic condensate is sensitive to the impurity scatter- 
ing @, 15|. Since the bound electron and hole carry the 
same momentum p, any scattering process that changes 
the momentum of electron and hole not identically, i.e., 
p — > p e for electron and p — > for hole, so that 
p e 7*= P/i, breaks the electron-hole pair (see Fig.[3|). This 
is the case for any impurities with the range of the scat- 
tering potential less than the interlayer distance d, since 
the potential of such impurities differs in the two layers. 
The effect of the impurity scattering on the excitonic con- 
densate was studied analytically for conventional systems 
in Refs. 0, EH and the theory is analogous to Abrikosov- 
Gorkov's theory for magnetic impurities in superconduc- 



tors. This approach has been very recently applied to 
graphene in Ref. [l6| ■ The main result of this study is 
that sufficiently short-range impurities with the scatter- 
ing time r destroy the excitonic condensate completely 
as soon as 



K/t > T c , 



(13) 



where T c is the transition temperature of the ideally 
clean system. Equivalently, for the condensate to exist, 
electron momentum has to be conserved at the scale of 
the correlation length hv/T c . Since the mean free path 
vt ~ 1/im of the order of the typical size of graphene 
samples corresponds to h/r ~ IK, considering the values 
T c max ~ lmK obtained, the condensate should be com- 
pletely suppressed at any temperature T < J 1 ™ 21 * even in 
the ballistic samples due to the boundary scattering. 

In conclusion, we have studied the possibility of the 
excitonic condensation in double-layer graphene systems. 
We have demonstrated that in order to properly deter- 
mine the transition temperature, it is essential to take 
the screening of the coupling interlayer Coulomb inter- 
action into account. The specifics of the graphene spec- 
trum (chirality and valley degrees of freedom) leads to 
a smaller interaction strength than in conventional semi- 
conductors and to an extremely small value < 1 mK of 
the transition temperature. This makes graphene-based 
systems disadvantageous for the observation of the exci- 
tonic condensation. 

After the present work was completed, we became 
aware of Ref. [17|, in which Eqs. (fTTj) and (fl"2")) were ob- 
tained. However, contrary to our main argument, that 
these equations are valid and provide good quantitative 
prediction also for moderate to strong Coulomb interac- 
tions (r s > 1), it was stated in Ref. [l7| that they should 
apply in the weak coupling limit (r s <C 1 or PFd ^> 1) 
only, whereas for r s ~ 1, analogously to Refs. [eOH, one 
could expect high T c of the order of room temperatures. 

We thank Anatoly F. Volkov and Sergey V. Syzranov 
for useful discussions. Financial support of SFB Transre- 
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